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Abstract 

O We investigate both bound and resonance states in 12 C embed- 



ded in a three-a-cluster continuum using a three-cluster microscopic 
model. The model relies on the Hyperspherical Harmonics basis to 
enumerate the channels describing the three-cluster discrete and con- 
tinuous spectrum states. It yields the most probable distribution of 
the three a-clusters in space, and the dominant decay modes of the 
three-cluster resonances. 

1 Introduction 

The 12 C nucleus is an interesting example of the so-called Borromean nuclei, 
. ^ as it has no bound states in any two-cluster subsystem of its three-cluster con- 

figuration. The lowest dissociation threshold (7.276 MeV above the ground 
03 state) is that of a three a particles disintegration. This three-cluster configu- 

ration is thus responsible to a great extent for the formation of a few bound, 
and many resonance states. The next threshold is of a two-cluster nature: 
ll B +p [lj. It opens when the excitation energy of l2 C exceeds 15.96 MeV. 
One therefore expects only a negligible influence of the latter channel on the 
bound and resonance states of 12 C in the vicinity of the a + a + a threshold. 

The l2 C nucleus is unique because of its excited "Hoyle state". This 
state is important in the context of the nucleosynthesis of carbon in helium- 
burning red giant stars. It is a + state with an energy of 7.65 MeV above 
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the ground state, or 0.4 MeV above the three-cluster a + a + a threshold. 
Its width is only 8.5 eV, indicating a long lifetime. One immediately relates 
this to the + state in 8 Be described by two a particles, with an energy of 
0.092 MeV above the a + a threshold, and a width of 5.57 eV. 

Many efforts have been made to reproduce the experimentally observed 
structure of 12 C, and to explore and understand the nature of the ground, 
excited and resonance states. This was e.g. done within so-called semi- 
microscopic models (considering structureless a-particles) [21 El HI El El [3 |8] 
and within fully microscopic models 13111011111113113111211211131131131121 

12011211 1221 123 123- 

A rather general feature of the calculations is that, with potentials which 
adequately reproduce the a-a interaction (this includes the phase shifts for 
+ , 2 + and 4 + states, and the position of the corresponding resonance states), 
one obtains a noticeably overbound ground state for 12 C. 

To determine the energies and widths of the resonance states created 
by a three-cluster continuum, only a few methods can be used. One popular 
method for obtaining the resonance properties in many-cluster, many channel 
systems is the Complex Scaling Method (see reviews [251 126] and references 
therein). Other methods start from a calculated form of the ^-matrix in a 
wide energy range, and determine the resonance states as the pole(s) of the 
^-matrix. The advantage of these methods is that they provide the scatter- 
ing quantities (such as phase-shifts, cross-sections, . . . ) and the resonance 
properties (energies and widths), as well as the wave functions of scattering 
and resonance states. The latter then allow one to obtain more information 
about the nature of the resonance states. 

12 C is known from theory and experiment (see, e.g., [32] and [33]) to 
have some very narrow resonances above the three a threshold. One may 
wonder why a system with several open channels does not decay instantly, 
but manifests these narrow resonance states. There are two possible answers 
to this question. First, a resonance state appears in one single channel of 
the multi-channel system. Such particular channel is usually weakly coupled 
to a number, or all, of the other open channels. It is well-known that this 
weak coupling of channels predetermines the existence of long-lived resonance 
states. Second, a resonance can be more or less uniformly distributed over all 
open channels, and the compound system needs (some) time for the resonance 
to be accumulated by one or a few number of open channels in order to decay 
into. Such a distribution over many open channels leads to very narrow 
resonances, as was predicted by A. Baz' [31]. It is referred to as diffusion-like 
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processes in scattering. This type of resonance is attributed to the effect that 
"the system spends most of its time wandering from one channel to another" 

In this paper we wish to calculate and analyze the bound and contin- 
uum structure of 12 C, and gain some insight in the nature of these states. 
Indeed, in some publications (e.g. (23 ESI ESI ED]) the suggestion for a dom- 
inant linear, chain-like, three-cluster structure appears for some of the 12 C 
resonances. We will look for confirmation of this structure. To this end, we 
determine the most probable configuration of the three a particles both in 
coordinate and momentum space. We also qualify those channels on which 
the resonance states of 12 C preferentially decay. 

The main results of this paper are obtained by applying the "Algebraic 
Model in a Hyperspherical Harmonics Basis" (AMHHB) [SI ESJ [36] on a con- 
figuration of three a-particles. In this model the three clusters are treated 
equally, and their relative motion descrybed by Hyperspherical Harmonics. 
The latter enumerate the channels of the three-cluster continuum and al- 
low to implement the correct boundary conditions for the three-cluster exit 
channels. The AMHHB has been applied successfully to study resonances 
in nuclei with a large excess of protons or neutrons such as 6 He, 6 Be, 5 H. 
The method provides the energies and widths of the resonances, and their 
total and partial widths, as well as the corresponding wave functions. The 
latter allow to analyze the nature of the resonance states. The results of this 
model are compared to those obtained in other, more or less comparable, 
microscopic descriptions from the literature, and to experiment. 

In the next section we elaborate on the method used to calculate the spec- 
trum of 12 C. Section three focuses on the results obtained using this method. 
We also present correlation functions and density functions to characterize 
more precisely the spatial configuration of the three a particles for specific 
resonance states. We also compare the results to those of other microscopic 
calculations as well as to experiment. 

2 The microscopic cluster model 

In this section we describe the microscopic model used to determine the 
structure of 12 C in the present paper. As it has already been introduced and 
used in several publications, we will limit ourselves to the most important 
notations and aspects of importance to the current calculations. 
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2.1 The three-cluster AMHHB model 



The three-cluster "Algebraic Model in a Hyperspherical Harmonics Basis" 
(AMHHB) [32 |35l E3 |36] will be applied to a single 12 C = a + a + a 
three-cluster configuration. 

This model takes a Hyperspherical Harmonics basis (HH) to characterize 
and enumerate the different three-cluster channels. In each of these channels 
an oscillator basis describes the radial behavior, and is used to expand the 
many-particle wave function. A matrix version of the Schrodinger equation is 
obtained after substitution of this wave function. It solved by the Algebraic 
Method (also called the Modified J-Matrix method [32]) for both bound and 
scattering states using the correct asymptotics. 

A similar approach, using the Hyperspherical Harmonics, was proposed 
in [39j HQ] in coordinate representation, using the generator coordinate tech- 
nique to solve the corresponding Schrodinger equation. 

The AMHHB wave function for 12 C is written as 

* = A{^(a l )^(a 2 )^(a 3 )f(ic,y)} (1) 
= A{^(a l )^(a 2 )^(a 3 )f(p,6;%y)} 

Yl C n P ,K,h,h\ n P^ K J^ l 2;LM;(p,e;%y)} 

n p ,K,h,l2 

where \n p , K, / 1; l 2 ; LM) is a cluster oscillator function [34] : 

\n p ,K,h,l 2 ;LM) = (2) 
A{$ (ax) $ (a 2 ) $ (a 3 ) R np>K (p) X K, h ,h W {Y h (x) Y h (y)} LM } 

These functions are enumerated by the number of hyperradial excitations n p , 
hyperspherical momentum K and two partial orbital momenta li, l 2 - The vec- 
tors x and y form a set of Jacobi coordinates, and p and 9 are hyperspherical 
coordinates, related to the Jacobi vectors by: 

p = a/x 2 + y 2 
|x| = pcosO, |y| = psin# (3) 

The notation x and y refers to unit length vectors. Vector x corresponds 
to the distance between two selected a particles, with an associated partial 
orbital angular momentum l 2 . Vector y is the displacement of the third a 
particle with respect to the center of mass of the other two, with an associated 
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angular momentum lx. The three quantum numbers c = {K, lx,h} determine 
the channels of the three-cluster system in the AMHHB. 

The fact that all three clusters are identical leads to some specific issues. 
The wave function Q for 12 C is antisymmetric with respect to the permu- 
tation of any pair of nucleons. Because the three clusters are identical, this 
function should be symmetric with respect to the permutation of any pair of 
alpha particles. This imposes constraints on the allowed quantum numbers 
of the wave function. Due to this symmetry, for instance, the partial orbital 
momentum l 2 of a two-cluster subsystem can only have even values. As the 
parity of 12 C states is defined as n = (— 1) 1l+1 ' 2 , it is fully determined by the 
partial orbital angular momentum l\ of the relative motion of the remaining 
cluster with respect to the two-cluster subsystem. 

In jH] and p2] it was suggested to use a symmetrization operator to 
construct the proper basis states. For a discussion on the symmetry of a 
system with three identical clusters we refer to [42] . 

The symmetrical Hyperspherical Harmonics basis for a three-particle sys- 
tem was realized many years ago (see, e.g., [43l E]). An explicit form of a 
few basis functions for small values of the total angular momentum (L = 0, 1 
and 2) can be derived. However it is extremely intricate to use for explicit 
calculation of matrix elements. 

An alternative approach to obtain such matrix elements without an ex- 
plicit realization of the basis functions consists in using the generating func- 
tion technique. One can indeed construct a generating function for the over- 
lap and hamiltonian kernels of 12 C, using the procedure explained in [34], 
that satisfies all required symmetry conditions, including the cluster sym- 
metric permutation behavior. Explicit matrix elements of the operators can 
then be obtained by using recurrence relations. The standard approach in 
the AMHHB is to extract matrix elements characterized by explicit 
quantum numbers. These, however, do not yet correspond to the desired 
symmetrical Harmonics. Indeed, the states \n p , K,li,l 2 ; LM) for fixed n p 
and K do not belong to the desired symmetrical irreducible representation 
of S(3), the permutation group of the three a clusters, with Young tableau 
[3]. They are, in fact, linear combinations of the Young tableau [3] and the 
non-symmetrical Young tableaus [2,1] and [111]. 

The antisymmetrization operator in the standard AMHHB basis has non- 
zero matrix elements 



( 



n p ,K,h,l 2 ;LM A n p ,K,h,l 2 ;LM 



) 



(4) 
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for fixed oscillator shells with N s h=2 n p + K = 2 n p + K. By selecting only 
the matrix elements with hyperradial quantumnumber n p = n p and hypermo- 
mentum K = K, one obtains relatively small matrices whose eigenfunctions 
K, u; LM) with non-zero eigenvalues are of the correct symmetrical Hy- 



\n 



p- 



perspherical type, due to the symmetry properties of the generating function. 
This procedure is similar to the procedure of obtaining the Pauli allowed 
states in three-cluster systems (for details see |45j). 
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Figure 1: Matrix elements of the antisymmetrization operator in the unsym- 
metrized Hyperspherical basis. 



This is demonstrated in Figure [T] where the diagonal matrix elements of 
the antisymmetrization operator between the original Hyperspherical Har- 
monics are displayed for total angular momentum L = 0, for all channels up 
to K = 8. One notices that the matrix elements 



n p ,K,li =l 2 ;L = Q 



A 



n p , K, lx = l 2 ;L = 



(5) 



do not tend to unity as one could expect, but to some fixed values. Analysis 
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shows that these asymptotic values of ([5]) correspond to the weights of the 
symmetrized Hyperspherical Harmonics with Young tableau [3], within the 
original Harmonic. 

The eigenvalues obtained after diagonalization however, which are ma- 
trix elemts of symmetrized Harmonics, do display the correct asymptotic 
behavior, i.e. all tend to unity, as can be seen in Figure [2j 




Figure 2: Matrix elements of the antisymmetrizator in the symmetrized Hy- 
perspherical basis. 



In Table [T] we display both the total number of (original) nonsymmetrized 
and of symmetrized channels for different values of the total orbital momen- 
tum. The symmetrization significantly reduces the number of channels com- 
patible with the maximal value of Hypermomentum K max . Only even values 
of the partial orbital momentum Z 2 are considered because of the symmetry 
rules for two-cluster subsystems. 

In [3l] and [35] the importance and meaning of the effective charges Z c ^ 
were defined in the context of the AMHHB, and explicitly discussed for the 



7 



Table 1: Number of channels for unsymmetrized and symmetrized Hyper- 
spherical Harmonics. 





0+ 


2+ 


4+ 


1- 


3~ 


K 


14 


14 


14 


13 


13 


N ch ({K,h,l 2 }) 


20 


44 


54 


28 


42 


N ch {{K,v}) 


8 


16 


19 


9 


14 



6 Be nucleus in three-cluster configuration A He + p + p. The effective charge 
determines the asymptotic form of the three-cluster potential originating 
from the Coulomb interaction, which has the form 

Vg> = Z -f (6) 

It was shown that it is of crucial importance for implementing the correct 
boundary conditions for the three-cluster continuum states. 

The symmetrization influences the behavior of the effective charges. In 
Table [2] we display the effective charges for the + state of 12 C, calculated 
in the original, nonsymmetrized, basis of the Hyperspherical Harmonics, for 
K-max — 8. One easily verifies that they coincide with those calculated in 



Table 2: Effective charges for the J w = + state of 12 C. 





(0,0,0) 


(4,0,0) 


(4,2,2) 


(6,0,0) 


(6,2,2) 


(8,0,0) 


(8,2,2) 


(8,4,4) 


(0,0,0) 


28.81 


2.47 


3.49 


2.74 


-2.74 


0.87 


0.00 


1.04 


(4,0,0) 


2.47 


32.157 


-1.13 


3.95 


-0.31 


4.67 


0.00 


1.95 


(4,2,2) 


3.49 


-1.13 


31.35 


1.72 


-4.30 


0.00 


0.66 


0.00 


(6,0,0) 


2.74 


3.95 


1.72 


33.48 


-2.51 


4.63 


0.00 


0.45 


(6,2,2) 


-2.74 


-0.31 


-4.30 


-2.51 


34.29 


0.00 


0.62 


0.00 


(8,0,0) 


0.87 


4.67 


0.00 


4.63 


0.00 


34.29 


0.00 


-2.38 


(8,2,2) 


0.00 


0.00 


0.66 


0.00 


0.62 


0.00 


33.08 


0.00 


(8,4,4) 


1.04 


1.95 


0.00 


0.45 


0.00 


-2.38 


0.00 


32.41 



In Table [3] we display the effective charges in the symmetrized basis. Only 
four channels remain after symmetrization. In particular no K = 2 channel 
remains, so we ommitted these also in Table [2] even though they have a 
non-zero contribution. 
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Table 3: Effective charges for the J n = + state of 12 C for symmetrized 
channels 



(K,v) 


(0,1) 


(4,1) 


(6,1) 


(8,1) 


(0,1) 


28.810 


4.277 


3.880 


1.139 


(4,1) 


4.277 


30.556 


5.217 


2.301 


(6,1) 


3.880 


5.217 


35.990 


1.457 


(8,1) 


1.139 


2.301 


1.457 


31.450 



It goes without saying that the asymptotic form of the effective three- 
cluster potential which originates from the nucleon-nucleon interaction 



is also influenced by the symmetrization. This asymptotic component is 
very important for obtaining the correct values of the S matrix. We do not 
dwell on its explicit form here, but apply a procedure similar to that for the 
effective charges. 



2.2 Phases, Eigenphases and Resonances 

After solving the system of linear equation of the AMHHB model, we obtain 
the wave functions of the continuous spectrum states, and the scattering 
S'-matrix. We consider two different representations of the S-matrix. 

In the first representation, the elements of the S-matrix are described 
through the phase shifts 5ij and the inelastic parameters r\if 

Sij = rjij exp (2iSij) (8) 

of which one usually only analyzes the diagonal matrix elements by display- 
ing the 5a and rju quantities. In the second representation the S'-matrix is 
reduced to diagonal form, leading to the so-called eigenphases, which now 
represent the elastic scattering of the many-channel system in independent 
(uncoupled) eigenchannels: 

||S|| = H^lp 1 ||D|| \\U\\ (9) 

Here \\U\\ is an orthogonal matrix, connecting both representations, and ||D|| 
is a diagonal matrix with nonzero elements 

D aa = exp (2i8 a ) (10) 
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defining the eigenphases 8 a . 

The phases shifts 5u, inelastic parameters rju and eigenphases 5 a then 
provide sufficiently detailed information about the channels that are involved 
in the production of resonance states. The eigenphases are used to extract 
the resonance positions and total widths in the traditional way 



d 2 5„ 



dE 2 



0. 



E — E r 



d8 a 
~d~E 



-1 



(11) 



whereas the orthogonal matrix \\U\\ leads to the partial decay widths of the 
resonance (for details see, e.g., [36J). 



2.3 Correlation functions and density distributions. 

As we pointed out, the AMHHB model allows to calculate the scattering 
properties, but also to obtain the obtain wave function at any energy, in 
particular at the resonance positions. The latter is of the utmost importance 
to analyze the nature of the system at these energies. 

Within the AMHHB model the solution is fully expressed by the ex- 
pansion coefficients {C npjC } and the 5-matrix. The expansion coefficients 
{Cn p , c } determine both the total three-cluster wave function of a compound 
system as well as the wave function of the relative motion of three clusters 
/ (x,y) (see eq. Q). 

The latter contains all information on the dynamic behavior of the three- 
cluster system for bound as well as continuum states. It is interesting to 
note that these coefficients are identical in both the representations of the 
wave function in coordinate and momentum space, because of the Fourier 
transform properties of the oscillator states. The wave function / (k, q) in 
momentum space has arguments that are directly related to the coordinate 
representation: k is the momentum of relative motion of two clusters, whereas 
q is the momentum of the third cluster with respect to the center of mass of 
the two-cluster subsystem. 

We obtain the density distribution in coordinate space as 

D(x,y) = D(p,9) = J |/(x,y)| 2 dxdy (12) 
and the corresponding correlation function as 

C(x,y)=C(p,6)=x 2 y 2 [ \f (x, y)| 2 dZ dy (13) 
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directly from the wave function of relative motion / (x, y). Both the density 
distribution and correlation function in momentum space are obtained in the 
same way using the wave function of relative motion in momentum space 
/(k,q). 

In a calculation with N c h open channels, one obtains N c h independent 
wave functions describing the elastic and inelastic processes in the many- 
channel system. It is quite impossible to analyze all of these wave functions 
when many channels are open. Some principles have to be set up on how 
to select the most important wave functions. In [36] we formulated some 
criteria for selecting the dominant wave function of a resonance. We will use 
the same criteria in this paper to select the "resonance wave functions" . 

3 Calculations and results 

In the present calculations for 12 C we consider for the nucleon-nucleon inter- 
action the Minnesota potential [46]. The oscillator basis is characterized by 
an oscillator length b = 1.2846 fm, to minimize the ground state energy of 
the a particle using the above potential. 

Parameter u of the Minnesota potential is taken to be u = 0.94 in order 
to reproduce the phase shifts for a + a scattering, and the + , 2 + and 4 + 
resonances in 8 Be. The same parameters were used by Arai [17J. 

The s Be = a + a two-cluster substructure is of key importance in the 
description of 12 C. We present a + a resonance properties in Table |4j The 

Table 4: Resonance properties for s Be obtained with different methods. 





AMOB 


Arai p2] 


J 77 


E, MeV 


T, keV 


E, MeV 


T, MeV 


0+ 


0.022 


6.30 10~ 10 


0.03 


<10" 6 


2+ 


2.93 


1.51 


2.9 


1.4 


4+ 


12.55 


5.01 


12.5 


4.8 



AMOB model takes a set of oscillator functions to describe the intercluster 
behavior, and the Algebraic Model to obtain the phase shifts for a + a 
scattering (see e.g. [4TJ ) . We include a comparison to the results of of Arai 
in his paper on 12 C [17] . where he uses the "analytical continuation of the S 
matrix to the complex plane" method with Complex Scaling to obtain the 
resonance characteristics. 
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These results form a first test of the consistency of the different expansion 
methods used, applied to the two-cluster subsystem. Although quite similar, 
one still notices that the resonance properties of the two-cluster a — a system 
have a slight dependence on the method used. 



3.1 The Potential and Coulomb interaction in AMHHB 

In Figures [3] and [4] the diagonal matrix elements of the nucleon-nucleon and 
Coulomb interactions within the AMHHB model are displayed, again for 
channels up to K = 8. One observes that the nucleon-nucleon interaction 
creates a deep potential well with a long tail in the Hyperspherical coordinate. 
This tail reflects the asymptotic form of the potential, indicated in ([7]). The 
matrix elements of the Coulomb interaction indicate the magnitude of the 
Coulomb barrier, which is the main factor for generating the resonance states 
in 12 C. 




Figure 3: Diagonal matrix elements of Vnn between symmetrized Hyper- 
spherical Harmonics. 
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Figure 4: Diagonal matrix elements of Vc between symmetrized Hyperspher- 
ical Harmonics. 

3.2 Phase shifts and eigenphases 

In Figure [5] we show results of the AMHHB calculations for the 2 + state 
in terms of the symmetrical Hyperspherical Harmonic channels through the 
(diagonal) phase shifts Sa and the inelastic parameters r/u . 

The scattering parameters are obtained from a calculation with maximal 
Hypermomentum K max = 14. One observes from Figure [5] that for small 
energies the channels are totally uncoupled f» 1). A first 2 + resonance 
appears at E = 2.731 MeV, and is mainly produced in the first channel with 
Hypermomentum K = 2, whereas a second resonance at energy E = 3.113 
MeV is dominated by Hypermomentum K = 4. The inelastic parameters for 
the first two channels have a pronounced minimum at the energy of the first 
resonance, and a shallow minimum at the second resonance energy. Also, 
the first resonance displays a "shadow resonance" behavior in the second 
channel. This is a typical behavior for resonances in a many-channel system 
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Figure 5: Diagonal phase shifts and inelastic parameters for the J 71 " = 2 
state. 



(see, for instance, the detailed analysis of two-channel resonances in 5 He in 
[48] ). The minimum in the inelastic parameters indicates that the compound 
system is being reconstructed at this energy, and transits from one channel 
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to another. 

In Figure [6] we display the corresponding eigenphase shifts 5 a for the first 
three eigenchannels. One observes now that both resonance states are mainly 
associated with the first eigenchannel, and that the second eigenchannel only 
contributes marginally 
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Figure 6: Eigenphase shifts for J n = 2 + for the first three eigenchannels. 



3.3 Convergence properties 

A convergence study of the energies (and widths) for bound and resonance 
states should indicate whether the Hilbert space is sufficiently large for sta- 
ble and reliable results. The AMHHB model space is characterized by two 
parameters: the maximal value of Hypermomentum K m3X , and the maximal 
value of the Hyperradial excitation n Pmax . Usually the choice is a compro- 
mise between the convergence of the results and the computational burden. 
A set of Hyperspherical Harmonics with K mgx = 14 for even parity states, 
and K max = 13 for odd parity states, seems sufficient and remains compu- 
tationally feasible. This choice accounts for a large number of three-cluster 
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configurations or, in other words, for a sufficient number of inherent (trian- 
gular) shapes for the three clusters. We refer to [38] for examples of most 
probable triangular shapes for the Hyperspherical Harmonics from K = to 
K = 10. 

A first convergence test considers the + , 2 + and 4 + bound states of 12 C, 
shown in Figure [7] as a function of K max . One observes that the deeply bound 
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Figure 7: Convergence of the bound states in AMHHB. 



states (J 71 = + , 2 + ) require significantly less Hyperspherical Harmonics for a 
converged energy than the shallow, or weakly bound, state with J n = 4 + . At 
least all Hyperspherical Harmonics with K mSuX > 6 are required to bind the 
latter state, whereas for J w = + one already obtains binding with a single 
Hyperspherical Harmonic with K = 0. Figure [7] further demonstrates that 
the above choice of K max amply leads to sufficient precision for the bound 
states. 

In Table [5] we turn to the energies and widths of the + and 2 + resonances 
obtained with increasing number of Hyperspherical Harmonics. One observes 
that sufficient convergence of the resonances occurs at K max — 12. It is 
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furthermore interesting to note that these resonances already appear with 
reasonable energy and width values when only the lowest channel (K = 
for the + , and K = 2 for the 2 + state) is considered. This is a remarkable 
result for 12 C, as e.g. for 6 Be it was impossible to generate a + resonance 
with a single K — channel (see [55]). 



Table 5: Energy (MeV) and width (keV) of the low- lying resonances in terms 

of Kxaax. 





K 





4 


6 


8 


10 


12 


14 


0+ 


E 


0.40 


0.75 


0.74 


0.72 


0.70 


0.68 


0.68 




r 


205.08 


13.40 


11.79 


7.10 


4.35 


2.71 


2.77 


0+ 


E 


1.15 


7.34 


6.09 


5.55 


5.54 


5.16 


5.14 




r 


510.16 


897.64 


422.50 


539.21 


586.08 


534.33 


523.46 


2+ 


E 




3.28 


2.89 


2.83 


2.78 


2.74 


2.73 




r 




30.19 


13.07 


11.85 


9.95 


8.84 


8.75 


2+ 


E 




3.50 


3.27 


3.22 


3.17 


3.14 


3.11 




r 




274.51 


351.57 


308.29 


280.23 


263.80 


246.78 



In all calculations we have considered states with hyperradial excitation 
up to n pmax =70, which covers a large range of intercluster distances, and 
reaches well into the asymptotic region. 



3.4 Partial widths. 

In Table [6] we display the energy, the total width (T) and the partial widths 
(Tj, i = 1,2,...) in the corresponding decay channels for the even parity 
resonances, and in Table [7] for the odd parity resonances. 

Table 6: Partial widths of the even parity resonances in 12 C. Energy in MeV, 
wi dths in keV. 



17 


0+ 


2+ 


2+ 


4+ 


E 


0.68 


2.78 


3.17 


5.60 


r 


2.79 


9.95 


280.24 


0.55 




K = 


2.79 


K = 2 


6.11 


K = 2 


13.46 


K = 4 


0.23 


r 2 


K = 4 





K = 4 


3.84 


K = 4 


278.89 


K = 6 


0.15 


r 3 


K = 6 





K = 6 


<10" 5 


K = 6 


<10" 5 


K = 8 


0.16 
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Table 7: Partial widths of the odd parity resonances in 12 C. Energy in MeV, 
widths in keV. 
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i- 


3~ 


E 


3.52 


0.67 


r 


0.21 


8.34 




K = 3 


0.206 


K = 3 


8.34 


r 2 


K = 5 


0.002 


K = 5 





r 3 


K = 7 


<10~ 5 


K = 7 






One observes that in most cases only one or two channels are responsible 
for the decay of the resonance states. The remaining channels contribute 
negligibly, and the corresponding partial width does not exceed 10 -5 keV. 
Only for the 4 + resonance a significant distribution over multiple channels is 
apparent. 

One should note that, although the resonances are created by only a 
few channels, the role of the other, very weakly coupled, channels is still 
important. This can be seen from Table [5] for the first + resonance: it is 
indeed generated mainly by the channel with minimal Hypermomentum K = 
0, but modified substantially with increasing number of Hypermomentum. 
The same applies to the other resonance states. 

3.5 Correlation functions and density distributions. 

In Figure [8] we show the correlation function for the 12 C ground state, and 
observe that this state displays a compact spatial configuration, as it is ex- 
pected for such a deeply bound state. The most probable shape of the three 
a-cluster system is an almost equilateral triangle with a distance between 
any two a-particles of approximately 3 fm. 

The correlation function for the first + resonance state on the other 
hand, shown in Figure [9j shows a more deformed system with two a particles 
relatively close to one another (about 3.5 fm) and the third alpha-particle 
further away (approximately 5 fm). So 12 C features a prolate triangle as a 
dominant configuration for this state. 

One also observes on Figure [9] a small maximum for the correlation func- 
tion corresponding to an almost linear configuration of three a particles, 
two of them being approximately 4 fm apart, and the third 0.2 fm away 
from their centre of mass. However, the weight of this linear configuration is 
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approximately 6 times less than the weight of the prolate triangular configu- 
ration. Our calculations therefore do not agree with other authors advancing 
a dominant linear structure j27J EEJ EHJ 130] . 



In Figure 10 we display the correlation function of the first resonance 
state in momentum space. One observes a huge maximum corresponding 
to relatively slowly moving a particles. A small maximum corresponding to 
faster moving alpha-particles is also present. 
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1 2 3 4 5 6 7 



r(a-a) (fm) 

Figure 9: Correlation function for the first + resonance state of 12 C in 
coordinate space. 

3.6 Comparison to the literature 

We now compare the AMHHB results to the existing literature. In Table [8] 
we display the AMHHB results to those of Arai [TT] and Pichler et al. [9], 
both obtained by the Complex Scaling Method (CSM). The latter authors 
[9] use a somewhat different value for the parameter u in the Minnesota 
potential, and a different oscillator length b; because of this, different results 
are obtained for the bound states. 

Comparison with the results of Arai [17] indicates that the AMHHB 
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0.1 0.2 0.3 0.4 0.5 

k(a-a) (Tm" 1 ) 

Figure 10: Correlation function for the first + resonance state of 12 C in 
momentum space. 



model leads to resonance states with higher energy and smaller widths than 
those obtained with the CSM. This can be attributed to the difference in the 
methods, and to the different Hilbert spaces. Formally the Hilbert space of 
basis functions used by Arai [T7| is quite close to the one considered in the 
AMHHB. Actually, in the present calculations the partial orbital momenta 
l\ and h are restricted by the condition 

L < l\ + I2 < A'max 

so that, for instance, for total orbital momentum L = 0, they run from 
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Table 8: Bound and resonance states of 12 C obtained with the AMHHB 
model, compared to CSM results from the literature. 



Method 


AMHHB 


CSM-Arai 


CSM-Pichler et al. 


Reference 


Present 


paper 


m\ 


i 






E, MeV 


r, keV 


E, MeV 


T, keV 


E, MeV 


T, keV 


0+ 


-11.372 




-11.37 




-10.43 








0.684 


2.71 


0.4 


< 1 


0.64 




14 




5.156 


534.00 


4.7 


1000 


5.43 




920 


2+ 


-8.931 




-8.93 




-7.63 








2.775 


9.95 


2.1 


800 


6.39 




1100 




3.170 


280.24 


4.9 


900 








4+ 


-3.208 
5.603 


7.82 


-3.21 
5.1 


2000 








\- 


3.516 


0.21 


3.4 


200 


3.71 


360 


3" 


0.672 


8.34 


0.6 


< 50 


1.16 




25 




4.348 


2.89 


7.1 


5400 


11.91 




1690 




5.433 


334.90 


9.6 


400 









h = h = to Zi = Z 2 = 6 with K max = 14. Arai on the other hand, restricted 
himself with li,l 2 < 4. In [35], [36], [37] we observed the tendency that the 
more Hyperspherical Harmonics (thus the more channels) are involved in 
the calculation, the smaller the resonance energy and width becomes. This 
tendency is again confirmed by the present AMHHB calculations. Thus some 
reduction of the width of the resonances, observed in our calculations with 
respect to Arai [17J, can be attributed to the larger number of channels in 
our model. 

Comparing the AMHHB results to the Complex Scaling Model calcula- 
tions of Pichler et al. [9], one observes that both yield close results for the 
first and second + resonance states. 

On the whole one can conclude that there is consistency in the results for 
resonance properties in all three microscopic models. 

3.7 Comparison to experiment 

In Table [9] we compare the theoretical AMHHB results for 12 C to available 
experimental data. 
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Table 9: Bound and resonance states of 12 C obtained with the AMHHB 
model, compared to experiment. 



Method 


AMHHB 


Experiment 


Reference 


Present 


paper 


m 


J n 


E, MeV 


r, keV 


E, MeV 


T, keV 


0+ 


-11.372 




-7.2746 






0.684 


2.71 


0.3796 ± 0.0002 


(8.5 ± 1.0) x HT 3 




5.156 


534.00 


3.0 ±0.3 


3000 ± 700 


2+ 


-8.931 




-2.8357 ±0.0003 






2.775 


9.95 


3.89 ±0.05 


430 ± 80 




3.170 


280.24 


8.17 ±0.04 


1500 ± 200 


4+ 


-3.208 










5.603 


7.82 


6.808 ±0.015 


258 ± 15 


1" 


3.516 


0.21 


3.569 ±0.016 


315 ± 25 


3^ 


0.672 
4.348 
5.433 


8.34 
2.89 
334.90 


2.366 ±0.005 


34 ±5 



One notices that the first + resonance state (the Hoyle state) appears in 
the current calculations as a narrow resonance with an energy of 0.684 MeV 
and width 2.7 keV, which is considerably wider than the experimental Hoyle 
state (about 8.5 10~ 3 keV). This contrasts with the generally observed fea- 
ture of the AMHHB calculations that the calculated widths are significantly 
less than the corresponding experimental widths of the 12 C resonances. The 
discrepancies between the theoretical and experimental data have essentially 
two origins. The first one relates to the choice of the nucleon-nucleon in- 
teraction: it has been tuned to reproduce the phase shifts and resonance 
properties for alpha-alpha scattering. As a result it leads to overbound + 
and 2 + states in 12 C, and binds the 4 + state. The second one relates to the 
specific choice of three-cluster model and corresponding model space, as well 
as to the method by which the energy and width of the resonance states are 
obtained. 
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3.8 Optimizing the nucleon-nucleon potential 

In this paper we used a Minnesota nucleon-nucleon potential tuned to repro- 
duce the phase shifts for a — a scattering, as well as the 8 Be resonances. This 
however leads to overbound + and 2 + states, and a bound 4 + state. More- 
over, the obtained resonance structure for the 12 C three-cluster continuum 
deviates from the experimentally observed one, which can also be attributed 
to the specific choice of semi-realistic nucleon-nucleon potential. 

We therefore wish to discuss the dependence of the results to the choice 
of parameter u on the results. To do so we use different criteria to optimize 
this parameter. We first determine a value to reproduce the ground state 
energy of 12 C, followed by an attempt to reproduce the energy and width of 
the + Hoyle state. 

In Figure [TT] we display the ground state energy as a function of the 
parameter u, compared to experiment (dashed line). One observes that the 
ground state is reproduced with u = 0.910. One observes a monotonously 




0.91 0.92 0.93 0.94 0.95 



ll 

Figure 11: Energy of the ground state as a function of parameter u of the 
Minnesota potential. 
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decreasing linear dependence of the ground state energy on u within the 
selected range. For the Hoyle state position and width the dependency is 
less trivial, as is shown in Figure 12 One however observes that the value 



1.0 - 




0.91 0.92 0.93 0.94 0.95 

u 

Figure 12: Position and total width of the first + resonance state as a 
function of parameter u. 

u = 0.948 reproduces the position of the Hoyle state, and leads to a close 
match for its width too. 

The correlation functions for the ground state and Hoyle state obtained 
with their respective optimal values were very close to the ones obtained with 
the value u = 0.94 and displayed in Figures [8] and [9j so that the conclusions 
remain unaltered. 
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4 Conclusions 



In this paper we described the 12 C nucleus with a three-cluster microscopic 
model. 

The model correctly handles the three-cluster continuum, i.e. correctly 
implements the suitable boundary conditions, by using a Hyperspherical Har- 
monics basis. It leads to the scattering matrix S in many-channel space, and 
the energy, total and partial widths of the resonance states and their corre- 
sponding wave functions can be obtained. 

It was shown that the obtained resonances of 12 C agree well with other 
methods, and that the lowest resonances are generated by only a few num- 
ber of weakly coupled channels, leading narrow resonance states. The partial 
widths determine the most probable channels for resonance decay. Correla- 
tion functions and density distributions revealed the dominant shape of the 
three-cluster triangle configuration for the lowest bound and resonance states 
of l2 C . There were no indications of a prominent linear three-cluster struc- 
ture for the resonance states. 

It was also shown that it is impossible to fix a unique value for the u 
parameter of the Minnesota nucleon-nucleon potential to fit all desired phys- 
ical properties for 12 C, and for the disintegrating a particles. However the 
qualitative conclusions remained unaltered under slight adaptation of u. 

As a final conclusion we can state that the model is consistent with other 
microscopic models using the Complex Scaling methodology to determine 
the energy and total width of three-cluster resonance states. 
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